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Abstract. We develop direct and inverse scattering theory for one-dimensional 
Schrodinger operators with steplike potentials which are asymptotically close 
to different finite-gap potentials on different half-axes. We give a complete 
characterization of the scattering data, which allow unique solvability of the 
inverse scattering problem in the class of perturbations with finite second mo- 
ment. 



1. Introduction 

In this paper we consider direct and inverse scattering theory for one-dimensional 
Schrodinger operators with steplike finite-gap background, using the approach by 
Marchenko [26] , 

To set the stage, let 

d 2 

(1.1) H ± =-_ +ftt (a;), xeR, 

be two (in general different) one-dimensional Schrodinger operators with real finite- 
gap potentials p±(x). Furthermore, let 

d 2 

(1.2) H = -— + q(x), xeR, 

be the "perturbed" operator with real potential q(x) € L\ oc such that 

/>±oo 

(1.3) ± \q(x) -p±(x)\(l +x 2 )dx < oo. 

Jo 

That is, the potential q{x) has different asymptotic behavior on different half-axes, 
and we will call it a steplike potential by analogy with the case of two different 
constant backgrounds. 

The scattering problem for the operators (|1.2[) - (|1.3[) is classical and arises in 
various physical applications, for example, when studying properties of the alloy 
of two different semi-infinite one dimensional crystals. We refer to the recent work 
[18] for a more detailed discussion of the history of such problems together with 
further references to the literature. In addition, this scattering problem is of course 
important for the solution of the Korteweg-de Vries (KdV) equation with initial 
data in these classes. 
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For a constant steplike background, that is, p±{x) = c± (c+ =/= c_ some con- 
stants), this problem was completely solved in [21 El [7] j including applications to 
the initial value problem for the KdV equation (01101): and the asymptotic behav- 
ior of the solution for large time ([22]). The case, when the potential vanishes on 
the right half axis (p+(x) = 0) and is asymptotically periodic finite band on the left 
half axis, was considered in [HI [301 EH [33] . The initial value problem for the KdV 
equation was also solved in this case ([12]). and long time asymptotics can be found 
in [TH [551 H3 (see also [5] for the Jacobi operator case) . The case of one-periodic 
background p- (x) = p + (x) was studied in [T^HTUHS] and [2H1 US] (for the Jacobi 
operators the same problem was considered in [8] with the extension to different 
background operators in the same isospectral class given in [10]). 

However, despite the fact that several special cases are well understood by now, 
only very little was known about the general situation considered here. In fact, 
the various mutual locations of the respective background spectra (cf. the example 
on page [7] below) and background Dirichlet eigenvalues, produce a multitude of 
different cases. To illustrate this, we mention that only a classification of all possible 
singularities of the transmission coefficient at the boundary of the spectra would 
require 16 case distinctions (this is for example the reason why we formulate its 
properties in terms of the Wronskian of the Jost solutions in Lemma I5!51 II below). 
Our goal here is to find a complete characterization of the scattering data for the 
operator H, that will allow us to solve the inverse scattering problem and to prove 
the uniqueness of the reconstructed potential in the class (jl.3|) with the second 
moment finite. In particular, we will do this without any restrictions on the mutual 
location of the respective spectra and without any restrictions on the location of 
the Dirichlet eigenvalues. In this respect, note that for example in (T5J [TH [TS] 
it is required that the Dirichlet eigenvalues do not coincide with the edges of the 
continuous spectrum. In fact, it was demonstrated in [8] (for the case of Jacobi 
operators) that these cases give rise to a different behaviour of the scattering data, 
which does not occur in the constant background case. Furthermore, the inclusion 
(and understanding) of this case is important for applications to the solution of the 
initial value problem of KdV equation, since these cases are unavoidable under the 
KdV flow. We refer to [9] for a detailed discussion (in case of the Toda lattice) and 
we will give a brief outline for the KdV equation in Section [6l 

On the other hand, we should also mention that there are two things which we 
do not address here: First of all, one could relax our decay assumption and replace 
the second moment in (jl.3|) by the first moment. Our approach is crafted in such a 
way that this causes no principal problems. To keep our presentation more readable 
we have decidedd not to include this case at this point. Secondly, one could allow 
general periodic potentials, that is, an infinite number of gaps. Again there are no 
serious impediments to treating this case. 

Finally, let us give a brief outline of the present paper. We start with some 
preliminary notations and list some standard facts of the spectral analysis for the 
background Hill operators in SectionO Then we study the properties of the scatter- 
ing matrix for steplike operator, paying particular attention to analytical properties 
of its entries at the edges of the continuous spectrum of operator H (Section [3j. In 
Section 0] we derive the Gel'fand-Levitan-Marchenko (GLM) equations and obtain 
complementary estimates on their kernels (see also Appendix [SJ. In this section we 
also formulate our main result, that characterize the scattering data (Theorem (4?3j. 
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Then we discuss the unique solvability of the GLM equations, that allows us to 
solve the inverse scattering problem. Section [5] is the most important section of the 
present paper. Here we discuss the scheme of the solution of the inverse scattering 
problem and prove the uniqueness of the reconstructed potential (Theorem !5-3[) . 
Our approach is modeled after the generalized Marchenko approach, developed in 
[26j . Our final Section [6] contains some applications to the KdV equation. 



2. The Weyl solutions of the background operators 

Let H± be two finite-gap one-dimensional Schrodinger operators associated with 
the potentials p±(^jH- Let s±(z, x), c±(z,x) be sin- and cos-type solutions of the 
equation 



(2.1) f--^ + p±(x)J y(x) = zy(x), z e C, 
associated with the initial conditions 

(2.2) s±(z,0) = 4(«,0) = 0, c±(z,0) = 4(z,0) = l, 

where the prime denotes the derivative with respect to x. 

It is well-known that finite-gap Schrodinger operators are associated with the 
Riemann surface of a square root of the type 



(2.3) 



\ 



2r± 

J[(z-Ef), £± <E± <...<Ef r± , 



where r± £ N. Moreover, H± are uniquely determined by fixing a Dirichlet divisor 

X^=i(Mi i a j )) where ^ e \^2j-\ > ^2j\ an d a f e { — 1, 1}- We refer the interested 
reader to [Tl ll7l[TOll25j for relevant background information. The reader not familiar 
with this theory can always think of the special case of periodic finite gap operators. 
The spectra of H± consist of r± + 1 bands 

(2.4) a ± := [E$ , Ef] U • • • U [E%_^E%_ X \ U • • ■ U [E± ± , oo). 

Note that in the special case where p± is periodic, we have merged all colliding 
bands. Let 

M r± := { M ±,..., M ± ± } 
be the set of Dirichlet eigenvalues and set 

n;iiO*-M±) 



(2.5) g±(z) 



where the branch of the square root is chosen such that we obtain a Herglotz- 
Nevanlinna function, 

(2.6) lm(g ± (z)) > for Im(z) > 0. 



Everywhere in this paper the sub or super index "+" (resp. "— ") refers to the background on 
the right (resp. left) half-axis. 
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Let us cut the complex plane along the spectrum <j± and denote the upper and 
lower sides of the cuts by a± and tj± . The corresponding points on these cuts will 
be denoted by A u and A , respectively. In particular, this means 

/(A u ) :=lim/(A + ie), /(A 1 ) := lim/(A — ie), A S a±. 

ej.0 e|0 



Condition (12. 6|) then implies 

(2.7) y ff± (A u )=Im( 5± (A u ))>0 for A E a±. 

Next consider the Weyl solutions tp±(z,x) and ■t/j±(z,x) of (|2.ip which are de- 
termined up to a multiplication constant, depending on z, by the requirement 

^ ± (v)ei 2 (l ± ), 

2.8 

resp. i>±(z, •) E L 2 {R T ) 



for z E C\ct±. We will normalize them according to ip±(z, 0) = $±(z,0) = 1 such 
that 

ip±(z, x) = c±(z, x) + m±(z)s±(z 7 x), 

(2.9) 

resp. ip±{z, x) = c±(z, x) + m±(z)s±(z, x), 

where 

(2.10) Tl±(z) = - — 7 — — , m ± (z) = — , 

are the Weyl m-functions. In the case of periodic operators, ij)±{z, x) and tp±(z, x) 
are of course just the Floquet solutions. They are equal to the branches on the 
upper/lower sheet of the Baker- Akhiezer function of H±. 

It is well-known (see, for example, |25j). that m±(z) — m±{z) = z fg±(z) ~ 1 . 
Equations (|2.2p and (|2.10|) then imply that the Wronskian of the Weyl solutions is 
equal to 

(2.H) W{4>±{z)^±{z))=T9±{z)- X . 

where W(f,g)(x) — f(x)g'(x) — f'(x)g(x) denotes the usual Wronski determinant. 
The set of band edges is given by 

(2.12) da± = {E±,Et,...,E± ± }. 

For every Dirichlct eigenvalue ^ the Weyl functions might have poles. If /if is in 
the interior of its gap, precisely one Weyl function m± or rh± will have a simple 
pole. Otherwise, if fj,j sits at an edge, both will have a square root singularity. 
Hence we divide the set of poles accordingly: 

M± = {[J.f | fif E (E 2 j-i, E 2 j) and m± has a simple pole}, 

A/± = | ^ E (E 2 j-i, E 2 j) and m± has a simple pole}, 

M± = {nf\»fe{E 2j _ 1 ,E 2j }}. 

Clearly M r± = M± U M± U M± . Then we have 

m ± (z) = -^- (l + o(l)), rh ± (z) = 0(l), 
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for z — * p £ M±, 



m±(z) = 0(l), m±{z) = J^-(l+ (l)). 



for z — > /i € M±, and 

m ± (z) = (1 + o(l)) , m±(z) = -^== (1 + o(l)) , 

Vz — E v z — E 

for z — »• _E e M±. Here C± ,C± denote some nonzero constants. 

In particular, we obtain the following properties of the Weyl solutions (see, e.g., 

[UUHl^EIIlIM]): 

Lemma 2.1. The Weyl solutions have the following properties: 

(i) The function ip±(z 1 x) (resp. i)±{z,x)) is holomorphic as a function of z in 
the domain C \ (a± U M±) (resp. C \ (<r± U M±)), takes real values on the 
set R \ a±, has simple poles at the points of the set M± (resp., M±). It is 
continuous up to the boundary a± U a\. except at the points from M± and 

(2.13) ^±(A u ) = V^ ± (A I )=^(A T y, \Eo-±. 
For E € M± the Weyl solutions satisfy 

^±{z,x)=o( K -j2== S j , $±(z,x) = O ( J— J , asz^EeM±. 

The same is true for rp'±(z, x) and ip'±(z,x). 

(ii) At the edges of the spectrum these functions possess the properties 

V>± (z, x) — ip± (z, x) = O z — E^j near E^da±\M±, 

and 

ip±(z,x) + $±(z,x) = 0(1) near E £ M± . 

(iii) When z — > oo the following asymptotic behavior holdE: 

i/>±(z,x) = e ±iV * x (l + 0(z- 1/2 )) and $± (z, x) = e^^ x (l + 0(z _1/2 )) . 

(iv) The functions ^p± (A, x) form a complete orthogonal system on the spectrum 
with respect to the weight 

(2.14) dp±(A) = -L 5± (A)dA, 

Z7T1 

namely 



(2.15) f ip ± (\,y)i)±(A,x)d P± (\)=8(x-y), 

where 5{x) is the Dirac delta distribution. Here we have used the notation 

(2.16) <f f(X)dp ± (X) := f f(X)d P± (X)- I f{\)d P± {\). 

J (7± J (T± J CTj_ 



2 Here lm(y/z) >0as2GC\R + . 
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3. The direct scattering problem 
Consider the equation 

(3.1) (--^+q(x)jy(x)=zy(x), z S C, 

with a potential q(x), satisfying condition (jl.3|) . This equation has two solutions 
<fi±(z,x), the Jost solutions, that are asymptotically close as x — * ±oo to the 
Weyl solutions of the background equations (12. 1|) and can be represented as (see 

[nmni]): 

/•±oo 

(3.2) (j)± (z, a?) = ip± (z, a?) ± / K± (a?, y)^± (z, y) dy, 

where K±(x,y) are real- valued, continuously difFerentiable with respect to both 
parameters, and satisfy the estimate 

/±oo 
\q(t)-p±(t)\dt. 
•+y 

Here C±{x) are continuous positive functions, monotonically decreasing (and, there- 
fore, bounded) as x — > ±oo (see Appendix IX)) . For A S o± U cr^ a second pair of 
solutions of (|3.1[) is given by 

/"±oo 

(3.4) <f>±(\,x) = $±(A,z)± / K±(x,y)4>±(\,y)dy, Ae^U^. 



Note ip±(\,x) = xjj±(X,x) for A 6 er±. 

We see that, by formulas (|3~2]1 . l|3~3]) , ([371)) . and (|2~TT|) . 



(3.5) W(0±(A),0 ± (A)) =±g±(\y 



Unlike the Jost solutions, the solutions (|3.4|) exist only on the upper and lower cuts 
of the spectrum of the corresponding background, and cannot be continued to the 
complex plane. 

The Jost solutions <j>± are holomorphic in the domains C\ (a± U M±) and inherit 
almost all the properties of their background counterparts, listed in Lcmma l2.ll (i)- 
(ii). In order to remove these singularities we introduce 



S ± (z) 


■■= n (* 






vfeM ± 




S ± (z) 


■■= n <* 


-^) n / 










■■= n (* 


-4) n / 









(3.6) 5±(z):= |] (*-/,*) {{ \/z-^f, 



where Y[ — 1 if there are no multipliers, and set 

(3.7) 4>±{z,z) = S±(z)ct>±(z,x), <f>±{z,x) = S±(z)ct)±(z,x). 

Lemma 3.1. The Jost solutions 4>±{z,x) have the following properties. 
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(i) For all x, the function (f>±(z,x) considered as function of z, is holomorphic 
in the domain C \ (cr± U M±), takes real values on the set R \ o~±, and has 
simple poles at the points of the set M± . Moreover, (j)± is continuous up to 
the boundary o~± U cr± . 

(ii) 0±(A, x) - <j>±{\,x) = O (a/A - E) for E e da± \ M±, and 
0±(A,x) +(f>±(X,x) = 0(1) forEeM±. 

Proof. Proof of this Lemma follows directly from (O, (O, (HHJ), Lemma [27T1 □ 

Now introduce the sets 

(3.8) cr (2) := cr+ n c4 1} = clos(cr± \ cr (2) ), cr:=(7 + Uer_, 

where cr is the (absolutely) continuous spectrum of H and <r+ U erf 1 '' , respectively 
cr*- 2 ) are the parts which are of multiplicity one, respectively two. In addition to 
the continuous part, H has a finite number of eigenvalues situated in the gaps, 
(j t [Ct\(T (see, e.g., [32]). We will use the notation int(cr±) for the interior of the 
spectrum, that is, int(cr±) := o~± \ dcr±. 

The following example illustrates the various possible locations of the spectra 
together with the Dirichlet eigenvalues. 

Example 3.2. Let H + be the two-band quasi-periodic operator with the spectrum 
on the set er+ = [Ei, E2] U [E4, +00) and H~ be the three band operator with the 
spectrum er_ = [Ei,E 2 ] U [Ea,Ei] U [E 5 ,+oo), where Ei < E 2 < ■ ■ ■ < E 5 (cf. 
Figure [Lj. Then a = [Ex,E 2 ] U [£3, +00), a™ = [E 4 ,E 5 ], a { } ] = [E 3 ,E A ], and 
cr( 2 ) = [Ei, E2] U [i?5, +00). Denote by the Dirichlet eigenvalue for the operator 
ii_ , that belongs to the closed gap [E 2 , £3] and let /i^ be the Dirichlet eigenvalue 
of H + from the gap [E 2 , E4] . 




i' 2 



Er, 



Figure 1. Typical mutual locations of ct_ and cr_|_. 



We may encounter various mutual locations of these eigenvalues. For example: 

• f-f t 1 ! and ^tilh € {E 2 ,E 3 ) (i.e., ^ £ M±), 

• l4 = Ih e {E 2 , E 3 ), 

• /if = E 2 ,iq £ E 2 (i.e. fi+ e M+), 

• ill € [E 3 , E4} (the Dirichlet eigenvalue is situated on the spectrum of multiplicity 
one) and \i{ 7^ E 3 , 

• etc. 

All such mutual locations of the Dirichlet eigenvalues imply different properties of 
the scattering data and have to be studied separately. 
Let 

(3.9) W(z):=W(<t>-(z,-),<p+(z,-)) 

be the Wronskian of two Jost solutions. This function is meromorphic in the domain 
C \ a with possible poles at the points M + U M_ U (M+ n M_ ) and with possible 
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square root singularities at the points M+ U M_ \ (M+ n M_). Set 

(3.10) W(z) = W$-{z),$+{z)), W{z) = W{4>-{z),4>+(z)). 

Then the function W(A) is holomorphic in the domain C \ R and continuous up 
to the boundary. But unlike the functions W(z) and W(z), it may not take real 
values on the set R \ a and complex conjugated values on the different sides of the 
spectrum. That is why it is more convenient to characterize the spectral properties 
of the operator H by means of the function W, which can have singularities at the 
points of the set M+ U M_ . We will study the precise character of these singularities 
below. 

Since the discrete spectrum of our operator H is finite, we can write it as 

(i d = {Ai, ...,A P } cR\cr. 
For every eigenvalue we introduce the corresponding norming constants 

(3.11) (t^)~ 2 = J 4>±(Xk,x) dx. 

Note that outside the spectrum W(z) — vanishes precisely at the eigenvalues. 
However, it might also vanish inside the spectrum at points in da- U <9cr+ and we 
will call such points virtual levels of the operator H 

(3.12) a v :={Ee<r\ W(E) = 0}. 

We will show a v C daU {da^} (Ida^) in Lemma I3~3l below. All other points E of 
the set da + U da- correspond to the generic case W(E) ^ 0. 

Now we begin our study of the properties of the scattering matrix. Introduce 
the scattering relations 



u,i 



(3.13) T T (A)</>±(A,x) = (j> T {\,x) +i? T (A)0 T (A,x), X e a^\ 

where the transmission and reflection coefficients are defined as usual, 

(3-14) 

rp r\\ ^QMA),<MA)) W/(0 T (A),0 ± (A)) , 

T±(A) -W(^(A),0 ± (A))' R±iX) --Wfa{\), ± (A))' AGCT± - 
Their characteristic properties will be given in the following lemma. 

Lemma 3.3. For the entries of the scattering matrix the following properties are 
valid: 



I. (a) T±(A U ) = T±(X l ) and i?±(A u ) = ^(A 1 ) for X € o±. 

(b) M^=R ± (X) forXea^. 
T±(A) ± 

(c) 1 - |i? ± (A)| 2 = ^ |T ± (A)| 2 for X e a^ . 

(d) T ± (A) = 1 + 0{\X\~^ 2 ) and R±(\) = 0{\X\-^ 2 ) for A — > oo. 

(e) R±(X)T ± (X) +R T (X)T ± (X) = for X e cr (2) . 

II. The functions T±(X) can be extended analytically to the domain C\ (aUM± U 
M±) and satisfy 

(3-15) - ~\ = - .T 1 =: W^(z), 

7+(2).9+(z) T-{z)g-{z) 

where W(z) possesses the following properties: 
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(a) The function W(z) = S+(z)5-(z)W(z) is holomorphic in the domain 
C \ a, with simple zeros at the points Xk, where 

fdW ( 



(3.i6) {^ Xk n =^k%Y 

Besides, it satisfies 



(3.17) W(X U ) = W(X l ), Aecr and W(X) G R for X e R \ a. 

(b) The function W(z) = 5 + (z)5-(z)W(z) is continuous on the set C \ a up 
to the boundary a u U c 1 . It can have zeros on the set daU^da^ nda^) 
and does not vanish at the other points of the set a . If W(E) — as 
E S do-U (da^ C)do- W ), then W(z) = y/7^E(C(E) + o(l)), C(E) ^ 0. 
III. (a) The reflection coefficient R±(X), is a continuous function on the set 
int(a± l ). 

(b) If E € do~+nda- and E ^ o~ v then the function R± (A) is also continuous 
at E. Moreover 



(3.18) R±(E) 



-1 forE£M ±1 
1 for E e M± . 



Proof. The proof is based on formulas (|3.14p . (|3.5|) , (|3.2[) , and Lemma |3~T1 

I. The symmetry property (a) follows from formulas (|3.14|) , (|2.13|) , (|3.2|) . and 
(|3~4l) . To verify (b) observe, that 4> T (X, x) e R as A <E int^^). Together with 
(|3~14ll and j3J]) this implies (b). Now let A e int(er( 2 )). Then by (f3TT3]) 

\T±\ 2 W(<P T ^) = (\R ± \ 2 - l)W{c^ ± ,JZ) 

and property (c) follows from (|3.5[) and (|2.14[) . To prove (d) we use p.5[) . (|2.7p 
and (|3.2p . Then property (hi) of Lemma l2. II allows us to proceed as in the proof of 
[26| Lemma 3.5.1]to obtain the necessary asymptotics. The consistency condition 
(e) and the identity (|3.15[) . considered on int(cr^ 2 ^), can be derived directly from 
the definition (|3~14| . 

II. (a). Except for (|3.16[) everything follows from the corresponding properties 
of 4>±(z, x) and it remains to show (|3.16jl . If VT(Ao) = for some Ao G C \ a, then 

(3.19) 4>±(X ,x) = c ± ct> T {X ,x) 

for some constants c± (depending on Ao) and satisfying c_ c + = 1. In particular, 
each zero of W (or W) outside the continuous spectrum, is a point of the discrete 
spectrum of H and vice versa. 

Let 7± be the norming constants defined in (|3.11[) for some point of the discrete 
spectrum Ao- Proceeding as in [26] one obtains 

(3.20) w^±(Ao,0),— ± (A o ,O)j =J cb 2 ± (X ,x) dx. 
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Equalities (f3TT9|) and (f3T20|) imply 



1± 2 = T<4 / <f>? F (\ ,x)dx ± / 4>±(\ ,x)dx 
Jo Jo 

= T c\W (<MA , 0), ^<M*o, 0)) ± T4^±(A , 0), J^±(Aq, 0) 
(3.21) = c ± -^W(4>-(X ),4> + (X )) 

and, since c_c+ = 1, we obtain (|3.16p . 
Item (b) will be shown in Lemma [B. 41 

III, (a) follows from the corresponding properties of 4>±(z, x) and from II, (b). 
To show III, (b) we use that by (|3.14[) the reflection coefficient has the represen- 
tation 



(3 - 22) ^ (A) " _ ^(0 ± (A),0 T (A)) " ± W(X) 

and is continuous on both sides of the set int(cr±) \ (Mzp U M T ). Moreover, 



\R±W\ 



W / (^±(A),0 T (A)) 



W(X) 



where the denominator does not vanish on the set a±\a v . Hence R±(X) is contin- 
uous on this set since both the numerator and denominator are. 

Next, let E e d<r± \ a v (in particular W(E) ^ 0). Then, if E £ M±, we use 
(|3~22| in the form 



(3.23) R ± (X } = -1 =p 

1 ' U W(X) 



which shows R±(X) — > —1 since (p±(X) — <f>±{X) — > by Lemma [Ql fiiV This settles 
the first case in (|3.18[) . Similarly, if E £ M±, we use (|3.22p in the form 



(3.24) S ± WW^ ± (X) + M^MX)). 

W(X) 



which shows R±(X) — > 1 since S±(X) — > and </>± (A) + <fi± (A) = O(l) by Lemma I5TT1 
(ii). This settles the second case in (13.18[) as well. □ 

We note that the behavior of T±(z) near the boundary points of the spectra can 
be read off from 



(3.25, W = ^ = -iMMS^. 

1 ' ±W g±(z)W(z) S ± (z) W{z) 



4. The Gel'fand-Levitan-Marchenko equation 

The aim of this section is to derive the inverse scattering problem equation (the 
Gel'fand-Levitan-Marchenko equation) and to discuss some additional properties 
of the scattering data, that are consequences of this equation. 
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Lemma 4.1. The inverse scattering problem (the GLM) equation has the form 

/>±oo 

(4.1) K ± (x,y)+F ± {x,y)± K ± (x,t)F ± {t,y)dt = 0, ±y > ±x 
where 

(4.2) F± (x, y)= <f R ± (A) ^± (A, x)^± (A, y)dp ± (A) + 

+ [ |T T (A)|V ± (A,x)^ ± (A,y)dp T (A) 

P 
k=l 

Proof. Consider the function 

(4.3) G± (z, x, y) = T± (z)0 T (z, x)-0± (2, 2/).9± (^) - i>±(z, x)i/j± (z, y)g± (z) 

:=G±(z,x,y)+G±(z,x,y), ±y > ±x, 

where x,y are considered as parameters. As a function of z it is meromorphic in 
the domain C \ a with simple poles at the points Afe of the discrete spectrum. It 
is continuous up to the boundary a u U cr, except for the points of the edges of 
background spectra, where 

(4.4) G±{z,x,y) = 0((z ~ Ey 1/2 ) as Eeda+Uda^. 
Since for z — > oo we have 

T (z,x) = e^ x (l + 0(z- 1/2 )), 4> T (z,x) = e Ti ^ x (l + 0(z- 1/2 )), 
tf>±(z, y) = e^^l + 0(z^ 2 )), T ± (z) = 1 + O^' 2 ) 
and g±{z) = ^= + 0(z^ 1 ) then 

(4.5) G ± (z,x,y)=e ±i ^ ( - y -^0(z- 1 ), ±y > ±x. 

Consider a closed contour T e consisting of a large circular arc together with some 
parts wrapping around the spectrum a inside this arc at a small distance from the 
spectrum. By the Cauchy theorem 

G±(z,x,y)dz = ^2 Rcs A fc G±(z,x,y). 

Estimate (|4.5|) allows us to apply Jordan's lemma. Since by (|4.4[) the limit value of 
G± as e — > is integrable on a, and the function G± has no poles at the points of 
the discrete spectrum, we arrive at 

(4.6) / G ± (X,x,y)dX = V Res Afc G' ± (X,x,y), ±y > ±x. 

Moreover, the function G± also does not contribute to the left part of (|4.6p since 
Gl(X u ,x,y) = G'l(X\x,y) for A € and, therefore, § aW G'l(X,x,y)dX = 0. In 
addition, § a G± (A, x, y)dX = for x ^ y by (j2~T5"l) . 
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Next we study the contribution of the function G'±. Properties I, (b) and (c) 
of Lemma 13.31 imply that 

(4.7) |i?±(A)|<l for Xemt(a^), \R ± (X)\ = 1 for A € <r±\ 

Now we consecutively use (|3~T5|) . (|3~T5)) . (j3~14l . (|TT5)) . (|3~2j) . (03) and again (|TT5)) . 
to obtain 

' ; G' ± (X,x,y)d\= £ T ± (\)<p T (\,x)iP ± (\,y)d P± (\) 



^± (A)0± (A, as) + 0± (A, aOj V± (A, y)dp± (A) 
R± (A)^± (A, a;)V± (A, y)d P± (A) + & 4>± (A, a;)^± (A, 2 /)d j o± (A) 

± dtK±(x,t)lf R ± (X)tp± (A, t)V»± (A, 2 /)d j0± (A) + *(t - y) 

/•±oo 

(4.8) = F r ,±(x,y)± K ± (x,t)F rt± (t,y)dt + K ± (x,y), 

J X 

where 

(4.9) F r ,±(^y)={ fl±W±(A,aOV±(A,i/)<MA). 



On the set u^ - both the numerator and denominator of the function G± have poles 
(resp., square root singularities) at points of the set <t|/ n (M± U (e?o+ n da_)) 
(resp. cr^ F 1) n(Af T \(M T nAf ± ))), but multiplying them, if necessary, by the functions 



(|3.6p . we can avoid singularities. Hence, without loss of generality, we can suppose 
n (M r+ U M r _) = 0. Then, since ip±(\,x) € M as A e cr^, 



(4.10) ±{ G' ± (X,x,y)dX = ± [ ^ ± (X,y) %g#-%g#WA. 



27rijr tr w ±v 2m y^D.u r = v \^ w-(A) w(A) 

Property I, (b) of Lemma |3~31 and (|3.13p imply 

^tOVe) = T T (A)0± (A, x) - J^U T (A,x). 

Therefore, 



<f> T (X,x) 4> T (X,x) _ / 1 T T (A) \ T T (A)0±(A,x 



W(A) T^(A) V 1 ^^) r T(^W/ W(A) 

(4 11) , 2Re(^(A)W^)r T (A) r T (A) 0± (A,,) 

(4 ' n) ~ MX ' x) WW W(T) ' 

But by (|XT5l 

7^ X (A)W(A) = |T^(A)| 2 3t (A) € «, for A € , 
thus, the first summand in (|4. 1 1|) vanishes. And using W = (T^g^z)^ 1 we arrive at 
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Combining (1412]) . (|2T3]) . (1410]) . (l3~2"l) and gl]) we have 



(4.13) — (f> G ± (X,x,y)dX = F Ci± (x,y) + K ± (x,y) ± 



±00 



K ± (x,t)F c< ±{t,y)dt, 



where 
(4.14) 

(4.15) 



F c,±(x,y) '•= F ri± {x,y) +F h ±(x,y), 



Fh,±(x,y) := 



(l),u 
=F 



^±(X,x)^ ± (X,y)\T T (X)\ 2 dp T (X). 



To derive the part of the GLM equation kernel, that correspond to the discrete 
spectrum (for function F Ct ± index c means the part, corresponding to the continuous 
spectrum), we apply $TT}, (|3~TU|) . ([3~T!)]) . (|3T2"Tj) and to the right hand side of 
(H"6l). Then, 



Res Afc G' ± (X,x,y) = - Res 
^±(A fc ,a:)^ ± (A fe ,y) 



T (A,a;)'0 ± (A,2/) 



E 



V^'(A fe )c ± , fc 

r>±oo 



T4^(A) 

= - E (lk) 2 ^>±(^,^ ± (X k ,y) 



(4.16) 

where 
(4.17) 



-F d ,±(x,y)T / K±{x,t)F d> ±(t,y)dt, 

J X 



Fd,±{x,y) := ^2 (l k t ) 2 ^±(Xk,x)ijj±(X k ,y). 

Combining ([4~6]) . (|4~13l) . and (f4TT6|) we finally obtain (j4~2|) . □ 
As is shown in Appendix \X[ the kernel F±(x, y) of the GLM equation satisfies 

Lemma 4.2. The kernel of the GLM equation possess the following properties: 

IV. The function F± (x, y) is continuously differentiable with respect to both vari- 
ables and there exists real-valued function q±(x), with 

±oo 

± I (1 + x 2 ) | q± (x) | dx < oo, /or all a £ K, 



sucft that 

(4.18) |F ± (a;,y)|<C±(a:)Q±(a: + y) ; 



(4.19) 



<C±(x) 



1± 



x + y 



+ Q±{x + y) 





/■±oo 




(4.20) 






where 






(4.21) 




Q±(a:) := 



(1 + x 2 ) dx < oo, 



|g±(*)|d*, 



and C± (x) > is a continuous function, which decreases monotonically as 
x — > ±00. 



14 



A. BOUTET DE MONVEL, I. EGOROVA, AND G. TESCHL 



In summary, we have obtained the following necessary conditions for the scat- 
tering data: 

Theorem 4.3 (necessary conditions for the scattering data). The scattering data 
S = {r+(X), T+(A), A e af; R_(A), T_(A), A G a u J; 

(4.22) Ai, .. . , Ap G R\ (<r+ U cr_), 7 ± . . . , 7 ± G R+} 

possess the properties I-III listed in Lemma \3.S\ . The functions F±(x,y), defined 
in (|4.2p . possess property IV from Lemma \4-2\ 

In fact, the conditions on the scattering data, given in this theorem are both 
necessary and sufficient for the solution of the scattering problem in the class 
(|1.1|) - (|1.3[) . The sufficiency of these conditions together with the algorithm for 
the solution of the inverse problem will be discussed in the next section. 

As a consequence of the GLM equation and its unique solvability (see Lemma [5~Tj) 
and also formula (|A.19j) we note 

Corollary 4.4. Suppose q(x) satisfies (|1.3p . then q(x) is uniquely determined by 
one of the sets of its "partial" scattering data S + or <S_, where 

S± = {r±(X), A G al- |T T (A)| 2 , A G af" 

(4.23) A 1 ,...,A p eR\(a+Ua_), 7l ± ,..., 7p ± GM+}. 

The question about the characterization of such sets (necessary and sufficient 
conditions) for potentials from the class (| 1 . 3[) is rather complicated and is still 
open. 



5. The inverse scattering problem 

Let H± be two one-dimensional finite-gap Schrodinger operators associated with 
potentials p± (x) as introduced in Section^ Let S be given scattering data with cor- 
responding kernels F±(x,y) satisfying the necessary conditions from Theorem 14.31 

We begin by showing that, given F±(x,y), the GLM equations (|4.ip can be 
solved for K±(x, y). 

Lemma 5.1. Under condition IV, (|4.18|) . the GLM equations (|4.1[) have unique 
real-valued solutions K±(x, ■) G L 1 (x,±oo) satisfying the estimates 

(5.1) \K ± {x,y)\<C ± (x)Q ± (x + y), ±y>±x. 

Here the functions Q±(x) are the same, as in (|4. 2 1[) . and C±(x) are functions of the 
same type, as in Lemma \4-.S\ (i.e. positive, continuous and decreasing as x — > ±oo). 

Moreover, under condition IV, (|4. 19[1 . the functions K±(x,y) are differ entiable 
with respect to each variable and satisfy the estimate (|A.22[) , where the functions 
q± (x) are the same as in (|4.19p and the functions C± (x) are of the same type as 
in (I5.ip . Besides, 



(5.2) ± / (1 + x 4 ) 



2^ 



-^-K ± (x,x) 

ax 



dx < oo, Va G 
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Proof. The solvability of (|4.1|) under condition (|4.18[) together with the estimate 
(|5.ip follows from considerations completely analogous to those ones used in the 
proof of Lemma IA.3I (see Remark IA.4|) . To prove uniqueness, first note that the 
GLM equations are generated by compact operators. Thus it is sufficient to prove, 
that the equation 

/■±oo 

(5.3) f{x)± F ± (x,y)f(y)dy = 



has only the trivial solution in the space L 1 (x,±oo). The proof is similar for the 
"+" and "— " cases, hence we give it only for the "+" case. Let f(y), y > x, be a 
nontrivial solution of (|5.3|) and set f(y) = for y < x. Since F + (x, y) is real- valued, 
we can assume /(y) is real-valued. Abbreviate by 

(5-4) /(A) = J i/; + (X,y)f(y)dy 

the generalized Fourier transform, generated by the spectral decomposition (|2.15D 
(cf. [M]). Recall that /(A) € Lj oc (a^ U a\) and /(A) = 0(1) as A -> +oo. 

Multiplying (|5.3[) by f(x), integrating over R, and applying (|5.4[) and (|4.2[) we 
have 



|/(A)|^ P+ (A) + 2Re / i?+(A)/(A)^ P+ (A) 
(5-5) +^ (i)u 7(A) 2 |T_(A)| 2 ^_(A) + ^( 7 +) 2 (Jj + (\ k ,y)f(y)dy^j = 0. 



The last two summands in (|5.5j) are nonnegative since /(A) € R for A € <j_ and 
■0+(Afe,x) €E R. Ignoring the last one and proceeding as in |26[ Lemma 3.5.3] we 
obtain 

(5.6) 2 / (1 - \R + (\)\)\f(X)\ 2 d P+ (X) + f /(A) 2 |T_(A)| 2 dp_(A) < 0. 



T< 2 )." 

Here we used that 



(l-|i? + (A)|)|/(A)| 2 ^ + (A)=0 



by condition I, (b). Now, since \R+(X)\ < 1, p+(X) > for A e int(a( 2 )) and 
P-(A) > for A G hit (erf 1 ''), we conclude that 

/(A) = for A e cr (2) U = cr_ . 

The function /(2:) can be defined by formula (|5.4p as a meromorphic function on 
C \ <t+. By our analysis it is even meromorphic on C \ (jj 1 and vanishes on cr_, 
thus f(z) is equal to zero and hence also f{x). 

The estimate (|5.2[) follows by literally repeating the proof of Lemma IA.3I □ 

Next, define two functions 

d 2 

(5.7) q±{x)=T-T^K±tx,x)+p ± [x), i£l 
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and note that estimate ()5.2j) implies 



(5.8) ± \q±(x) -p±(x)|(l + x 2 )dx < oo, a£i 



Lemma 5.2. The functions 4>±(z,x), defined by 



(5.9) <t>±(z,x) = ip±{z,x)± I K±(x,y)il>±(z,y)dy, 



■J X 



solve the equations 

(5.10) ^-^+<?±(a;)^ <t>±(z, x) = zc/>±(z,x), 

where q±{x) are defined by (|5.7|) . 
Proof. Consider the two operator^ 



On the corresponding half-axes the potentials q± (x) are asymptotically close to our 
background potentials p± (x) . Define the integral operators 



Under the assumption, that the kernel F±(x,y) of the GLM equation is twice 
continuously differentiable, we infer from (14. ip that the function K±(x,y) is also 
twice differentiable (see the proof of Lemma fA.SI) . Moreover one can prove, literally 
following [14] , that the identity 



is valid. This identity implies (|5.10|) . To obtain equality (|5.10|> without assumption 
of existence of the second derivatives, one can literally follow the proof of [26| 



The remaining problem is to show q+{x) = q~(x) under conditions II and III 
on the scattering data S. 

Theorem 5.3 (uniqueness of the reconstructed potential). Let the scattering data 
S, defined as in (]4.22[) . satisfy conditions I, (a)— (d), II, III, (a) and IV. Then 
each of the GLM equations (|4.1|) has a unique solution K±(x,y), satisfying the 
estimate (|5.2|) . The functions q±{x), defined by (|5.7p . satisfy (|5.8j) . 

Under additional conditions III, (b) and I, (e), these two functions coincide, 
q—(x) = q+(x) —: q(x), and the data S are the scattering data for the Schrodinger 
operator with potential q(x). 

Proof. To prove the uniqueness of the reconstructed potential we follow the method 
proposed in [26] . Recall that, according to [25l[34], we have the inversion formula for 




H±JC± = IC± H± 1 



Theorem 3.3.11. 



□ 



^We don't know H± is limit point at =poo yet, but this will not be used. 
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the generalized Fourier transform, generated by the spectral decomposition (|2.15j) 
and applied to the function /(A) e L\ oc {&± U a l ±), /(A) = 0(1), A — * +00: 

f(y) = <f f(X)^±(X,y)d P± (X), 



(5-11) /(A)= / f(y)i>±{X,y)dy 



Split the kernel of the GLM equation (|4.2p according to F±(x,y) — F rt ±(x,y) + 

F h ,±( x >y) + F d,±( x >y) ( cf - (EH), 6331, (SUZD)- 

We begin by considering the following part of the GLM equation 

/>±oo 

(5.12) G±(x,y):=F r ,±(x,y)± K±(x,t)F r ,±{t t y)dt, 

J X 

where K±(x,y) are the solutions of GLM equations. By condition I, (b)— (c) we 
have |.R±(A)| < 1 and properties (iii) and (i) of Lemma l2.1l implv. that we can take 
/(A) = R±(\)<ip±(\,x) in 1(531]) . Using (gSJ) we obtain 



(5.13) / F r , ± (x,y)ip±(X,y)dy = R ± {X)^ ± (X,x) 

JR 

and (|3.2p consequently implies 



.1 



(5.14) G±(x,y)*l;±(X,y)dy = R ± (X)<j> ± (X,x), X e <t u ± > . 

On the other hand, invoking the GLM equations we have for ±y > ±x, 
G± (x, y) =-K± (x, y) - F hi ± (x, y) - F dt ±(x, y) 



±00 



T / i u dp T (A)|T T (A)|V ± (A,y) / X±(x,t)^±(A,i)A 
T^(7^) 2 ^±(A fe ,y) / K-±(a: > t)^±(Afc ) t)dt- 



fe=l 

Again using (|3.2| this gives 

G±(a;,tf) = -if±(x,»)- / |J 1 T (A)|V±(A,2/)^±(A, a ;)^ T (A) 

P 

(5.15) -^(7 t ± ) 2 V^(A fc ,y)0 ± (A fc , a; ). 
fc=i 

Now we use this formula to evaluate 

G±{x,y)i) ± (X,y)dy =T / G±(x,y)ip ± (X,y)dy =p / K±{x,y)i)±{X,y)dy 



± / (l)i JIV(Or^±(C,a:)W^±(€,a:),V±(A,a:)) A 



±E(7, ± ) 2 ^ ± (A fc ,^) 



iy(Vi ± (A fc ,a;),V» ± (A,a;)) 

fc =i Afc _ A 

Here we have used 

/*±oo 

(C-A) / ^±(A J »)V±(e,v)dl/ = W(^±(A,a:),^±(e,a;)), 
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which follows from Green's formula and the fact, that ip±(£,y) — * 0, i/}'±(£,y) 
for £ (fi o~± as y — ► ±00. 

Combining this formula with (|5.14[) and using (|3.4[) we infer the relation 

(5.16) R ± (A) <t> ± (A, x) + <f>± (A, 5) = T± (A)0 T (A, x), A € crj 1 , 
where 

(5.17) T (A,x) ^jT^t ($±(\x)T J G ± {x,y)$ ± {\,y)dy 

|T T (£)| 2 <M£, x)W<$± K, x), ^ ± (A, x))^^ 



f 



±E(7, ± ) 2 ^±(A fe , a; ) 



fc=i 



iy(^ ± (A fc ,x)> ± (A,x)) 
At — A 



It turns out that, in spite of the fact that 6^{\x) is defined via the background 
solutions corresponding to the opposite half-axis R±, it shares a series of properties 
with 0zp(A,x). Namely, we prove 

Lemma 5.4. Let 0zp(z,x) be defined by formula f|5 . 1 7[) on the set <J± ■ 

(i) The function 0^{z,x) = 5^(z)0^(z, x) admits an analytical extension to the 
domain C \ o~. 

(ii) The function 9zp(z, x) is continuous up to the boundary er 11,1 except possibly at 
the points dcr + U do~ . Furthermore, 



(5.18) T (X u ,x) 



:(A', x) G R, for A G M \ ct 



1 (^(A 1 , x), /or A G int(o"=p). 

(iii) For large z the function 9z^(z,x) has the following asymptotic behavior 

T {z,x) = e^ x (l + O(z- 1 / 2 )), z^oo. 

(iv) W(9zp(z, x), <j>± (z, x)J = ±W(z), where W(z) is defined by formula p,15p . 

Remark 5.5. Note that we did not establish the connection between the function 
W(z) and the function W((j) + (z,x),<f>_(z, x)), which can depend on x, because 
<j) + and <p- are the solutions of Schrodinger equations corresponding to possibly 
different potentials q+ and gL. 

Proof. To show (i) we will show that 0^(z,x) has an analytic extension to C \ a. 
We will study each term in (|5 . 1 7[) separately. 
For the first one, 

(5-19) ^ : =W' 

it is immediate by (|3.25| that (^(z, x) = S T (z)( T (z, x) has the required property. 
This also covers the second term since G±(x, ■ ) G i 1 (R) is real- valued. 

Now we discuss the properties of the Cauchy-type integral in the representation 
(|5.17p . Multiplying it by T± (z), we represent the third summand in (|5.17|) as 



(5.20) H T (z,x) -T^J^hrfr&x) 
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where the integrand, due to (|3.10p . (|2.14|) . and (13.15| . has the representation 

(5.21) = 1^(01 l*±(fll j> ± (£, x )W(ii; ± {Z,x),( T {z,x)). 

By property II, (b) the function W(£) has no zeros in the interior of a^' u . Thus, 
for z cr^\ the function h T (z, .,x) is bounded in the interior of and the only 
possible singularities can arise at the boundary. We claim 

fO(VFTE) tor E£a v , m 

(5.22) M^,*) = { {orEea ^ Ee64\^E. 

This follows from = 0(V^~E) together with W(£) = 0(1) if E $ a v and 

W(0 = C(E)(VT=h(l + o(l)), C(E) by II, (b) if E G a v . 

So hzf: is integrable and the third summand of (|5.17|) also inherits the properties 
of <&{z,x). 

Finally, let us consider the last summand in (|5.17j) . It again inherits everything 
from C=f(^i x ) except for possible additional poles at the points A&. However, these 
are canceled by the fact that the function W(z) vanishes for z = Xk- 

(ii). Next we look at the boundary values. The only nontrivial term is of course 
the Cauchy-type integral (|5.20|) as z — > A G mt(a^). First of all observe that by 
(|2~TTj) and pTT5)) we have 

W$± (A, x),ip± (z,x)) 



T±(z) 



±(5 ± W)(X), 



where the function S±W is bounded and non zero for A G int(4 1} ) by II, (a). 
Hence the Plemelj formula applied to (|5.20|) gives 

4>±(X,x) f h T (\,£,x) . ( i) )U 



H^(X,x)=± , +v ^p- de. AGint(a^' u ) 

TV 2S ± (X)g T (X)W(X) e-A ' T 

where both terms are finite. Here denotes the principal value integral. 
Hence the boundary values away from da + U <9cr_ exist and we have 



(5.23) 6 T (X u ,x) = 6 T (X l ,x), Xecr + U<T-. 
Moreover, by property I, (b) we have 

(5.24) e T =T± 1 (R ± <p ± +^) = ^ + ^€M. for A G *£\ 

T± J-± 

from which 

(5.25) 6 T (X u ,x) = 0zf(X l ,x), Xea { i\ 

follows. Combining (|5T23"1) and ([05]) we have (|5T8|) . 

(iii). For \z\ — > oo due to properties (iii) of Lemma 1 2. II and I, (d) we have 

(5.26) ( T (z,x) = c^ x (l + Oiz- 1 ^ . 
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Since the last two terms in (|5.17|) are 0(z^ 1 ) we obtain 

9 T (z,x) =e^^ x U + J G ± (x,x T ty^dt + Oiz- 1 / 2 )^ 

which implies (iii) since G±(x,y) is differentiable with respect to y. 
(iv). From (|5.16|) (invoking (|3.15p ) we obtain 

(5.27) W (9 T (z,x),<f>±(z, x)) = ±W{z) 

for z G a± . Hence equality holds for all z G C by analytical continuation. □ 

Corollary 5.6. The function 9^{z, x) admits an analytical extension to the domain 
C\a T . 

Proof. Property (i) holds for z G C \ a. Relation (|5.18p implies that 6^ has no 
jump across z G mt(a^). To finish the proof we need to show that the possible 
remaining singularities at E G da^ n da are removable. This follows from (cf. 

(£25)) 

(5.28) Ct (z, x) = —j= ^ Z) 5 ± (z)i>±(z, x) 

which shows ( T (z,x) — 0((z — E)^ 1 / 2 ) and hence 9 T (z,x) = 0((z — E)^ 1 / 2 ) for 
E G a ( l ] n da. 

However, let us emphasize at this point that the behavior of 9±(z,x) at the 
remaining edges is a more subtle question to be discussed later. □ 

Eliminating <f>± from 



R± (A) <j>± (A, x) + <j>± (A, x) = 9 T (A, x) T± (A) 



4 R± (A) <t>± (A, x) + <t>± (A, x) = 9 T (A, x) T± (A) 
we obtain 

0± (A, x) (1 - |i?±(A)| 2 ) = 9 T (X,x)T ± (X) - R^X)9 T (X, x) T± (A). 
Next, using I, (c), II and the consistency condition I, (e) then shows 



T T (X)J> ± (X,x) = 9 T (X,x) R±i ^ ±{X) 9 T (X,x) 

= 9 T {X,x) + R T (X)9 T (X,x), AGa (2) . 

This equation together with (|5.16[) gives us a system from which we can eliminate 
the reflection coefficients R±. We obtain 

(5.29) 

T 1 ±(A)(0±(A)0 T (A) - 0±(A)MA)) = 0±(A)0±(A) - 0±(A)0 ± (A), A G a^ u \ 

Now introduce the function 

(5.30) G(z) := G(z, x) = * +(z ' x)4> ^ ^ +( "' x) °- ^ x) 

which is well defined in the domain z G C \ (a U ad U M + U M_). 
From (j!T2"9"|) and |3~T5)) we see, that 

(5.31) G(A) = -(^(A)4W-p)ft i (A))ji(A) 1 Ae^K 
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So we need to study the properties of G(z,x) as a function of z (regarding i as a 
fixed parameter). Our aim is to prove that G{z 1 x) = 0. This will follow from 

Lemma 5.7. The function G{z,x) possess the following properties. 

00 

(5.32) G(A U , x) = G{X\x) e R for A e R \ {da- U 9(7+ U cr rf ). 

(ii) It has removable singularities at the points da- U da + U ad, where a e i :— 
{Ai, . . . , A p }. 

Proof, (i). We can rewrite G(z,x) as 

4>+(z, x)<f>_ (z,x) — 0+{z, x)8_ (z, x) 



(5.33) 



G(z, x) 



W(z) 



where 9±(z, x) = S±(z)9±(z, x) as usual. The numerator is bounded near the points 
under consideration, and the denominator does not vanish there. Thus G(z, x) has 
no singularities at the points (M+ U M_) \ 

Furthermore, by Lemma T5.4[ II, (a), and Lemma T3. 11 (i) we know that G(z,x) 
has continuous limiting values on the sets er_ and er+ , except possibly at the edges, 
satisfying 

G{X u ,x) = G(X\x), Xea+Ua-. 

Hence, if we can show that these limits are real, they will be equal and G(z, x) will 
extend to a meromorphic function on C, that is, (i) holds. 

First of all observe that from 1(535]) . (jOTj) . (|2~7|) . and Lemma |3~T1 (i). it follows, 
that 

(5.34) G(A U , x) = G(A 1 , x) € M, A E int(cr^). 
Thus, it remains to prove 

(5.35) G(X u ,x) = G(X\x) e R for A G int^) U int^). 

So let us show that there is no jump on the set mt^a^) U int(oj^). Abbreviate 

(5.36) 



[G] := G(A) - G(A) = 



0+0- 




'9 +9 J 


w 




w 



A e 



and let us drop some dependencies until the end of this lemma for notational sim- 
plicity. 

Let A e int(o-£ ) ' u ), then 0±,0± £ R and T T = (W g^)" 1 . Thus, by condition 
(I), (b) and (|5TT6|) for A e int (a^ ] ) we obtain 



(5.37) 



0+0- 




[0+1 


w 


= 0± 





-#+0± 



= T T ) =-M±0±l I Vl 



Since g± € R for A € int(<T^ ),u ), jgjj) implies 



TL 



-5± [^2±] 



The only non-real summand in (|5.17[) is the Cauchy-type integral. The Plemelj 
formula applied to this integral gives 

[0 T T ± ] =±g T 0±\T T \ 2 W(4> ± J ± ) = 5 +0±|T+| 2 — ■ 

9± 
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Thus by ([537]) 
(5.38) 



'6+9-' 




(f>+<f>- 


W 




w 



-9t<I>±0±\T t \\ Aeint(aW). 



Since W ^ and s T ^ for A G int(a^'), the function 



,2 _ ^ 0±0± 



St \W\ 2 

is bounded on the set under consideration. Finally, (|5.38|) and (|5.36|) imply (|5.35|) . 

(ii) Now we prove, that the function G{z, x) has removable singularities at the 
points <9cr_ U da + U ad- Divide this set in four subsets 

(5.39) Qf = <9cr (2) n int(cr T ), = da { ± ] n da, fl 3 = da^ n da +1 and ft 4 = a d . 

In our example we have f2]*~ = 0, fi^ = {E 5 }, Qt — 0, = {-E3}, and f2 3 = 
{Ei,E2, E4}, 

By condition II, (b) all singularities of G are at most isolated poles. Thus it is 
sufficient to check that 



(5.40) 



G(z) = o((z-E)-i) 



from some direction in the complex plane. 

• f2i: Consider E G fi^ (the case E G fl^ is completely analogous). We will 
study lim\^E G(\,x) as A G int(er( 2 )) using identity (|5.3ip . We have 0_ = O(l), 
g_ = O(l), and W(E) 7^ 0. Moreover, from Lemma [37X1 respectively II we deduce 



<MA) = 
which shows 




E<£M+, 
EEM+, 



1 



o 



T+(X) \ (1), 



E4M. 



+) 



EeM+, 



? _ (A) = 0+(A) +J R + (A)^ + (A) = ^ 



1 



1 



and finishes the 



for A G a^ . Inserting this into (|5.31[) shows G(A, x) = O 
case Bgfli. 

• SI2 : Suppose now that E G 9<r_ Plc?<7 (the case E G 9o+ Plc?<7 can be treated in 
the same manner). Now we cannot use (|5.31[) . so we proceed directly from formula 



(|5.30p estimating the summands 
(b) we have 



w 



(5.41) 



W 



w 



and + separately. By Lemma |3. II and II, 



= O 



E 



Hence the first term is under control and it remains to investigate the second one. 
We investigate the limit of G from the set int(<7_ ). By (|5.30|) 



(5.42) 



II" 



_T_«?_ = (</)_+0_i?_)</_ 




E e M_, 
E £ M-, 



A G int(cr_). 



Next we will estimate 6>_ using 1|5.17|) . First, let E £ a v , that is W(E) 7^ 0. Using 
the same notation, see (|5.20[) . as in the proof of Lemma l5~4l we can split 0-(A) 
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according to 

(5.43) 0-(\) = x (A)+0 2 (A), 
where 

(5.44) 2 ( A ) = -L/ h _ M ^L, l (X) = 0.(X)-0 2 (X). 

2tti J a W- u £ - A 

We see that (cf. ([OS]) ) 

9l= o( ( _)4o (1 ,=H^)- 

(0(1), £^Af_, 
since E ^ a + . Next, we use (|5.44[) . where (cf. (|5.2ip ) 



h-(X,0 = ^ g C(OQ(C-(A)) 
|W(OI 

with C(£) some bounded function near E. Hence 02 = 0(£_) as well, which implies 
together with (jS"!^ that 

9+(\)0-(\) = Q ( 1 



W(X) \VX^E 



This finishes the case E £ <j v . 

Now let E E <r v . As before we have 



(5.45) 0i=O(C-] 



O(l), E e M_, 



(^(VA - ^), £^ M_ 
For the Cauchy-type integral we now have 

MA,0 = 7f=Lo(C-(A)) 
and [27, Eq. (29.8)] implies 



(5.46) 2 (A) 



.o(l), Ei : M 



Combining (|5.45[) . (|5.46[) . and (|5.42|) finishes the second case. 

• £13: The first step is similar as in the case of f?2- In particular, 1|5.41|) is valid 
for E € f2 3 and z € C close to E. Estimate (|5.42|) is valid for A € int(cr_). The only 
difference being that in estimate for has an additional square root singularity 
since E € 9cr+. That is, instead of (|5.45[) and (|5.46|) we now have 

M*) = C-(*)(c + °(i)) = 

However, this is not good enough unless we can show C = 0, in which case the big 
O will turn into a small o and we are done. It is sufficient to show C = from one 
direction, say A € o+, which can be done using the scattering relations for 9- as 
follows. 
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If E € a v this follows directly from 
g _ (A) _ ^ + fl+(A)<MA) _ Q ( W(X)$+(X) 




A^M + , 



T+(X) \5-{\)^\~E J 

Otherwise, if E (£ a v , then we have two representations 

(5.47) 0_(A) = ((07(A) - </>+ (A)) + 0+(A)(i? + (A) + 1)) 

(5.48) 0_(A) = _^y ((07(A) + + (A)) + + (A)(i? + (A)-l)), A G M+. 

For (|5T4Tj) we use ^+(A) - </>+(A) = o(l) by Lemma GO (ii) and i?+(A) + 1 = o(l) 
by condition III. Now one checks 



r+(A) 



+ 



0+(z) 



) n 7 \ - A £ M+, 



which shows G(z, at) = o( for A ^ M+. For (pT48f we use 4> + (\)+<p + (\) = 0(1) 
and i?+(A) — 1 = o(l) as A S M+. Since in this case T ^ z ) = O ( ^-b ) an< ^ 




3+(*) O -O, , EiM., 



this implies G(z,a;) = o(-rrg) for A e Af + as required. This finishes the case 

E e n 3 . 

• SI4: Finally we have to check, that the singularities at the points of the discrete 
spectrum are also removable. Since W(z) has simple zeros at points A&, it suffices 
to check that 

(5.49) 6+(\ k ,x)L{\ k ,x) = 4>-(\ k ,x)0 + {\ k ,x). 

Lemma [5T41 shows . that 9 T := S T T , defined by (|5.17[) . are continuous at the points 
M±. Since ( ^T± 1 )(A fc ) = and (<5 T T± 1 ^±)(A fe ) = 0, then the only the last 
summand in (|5. 1T[) is non-zero. Computing the limit of this summand at X k and 
using (|3TTD|) . ([3T4| we obtain 

(5.50) ^ (Afe ) = ^^l( 7± )^ ±( A fc ), 

which together with (|3 . implies (|5.49[> . □ 

Lemma 15.71 implies, that G{z, x) is an entire function for fixed x. Since in ad- 
dition G{z,x) — > as z — > co, Liouville's theorem implies G(z,x) = for every x. 
Therefore we have the equalities 

(5.51) (/> + (z,x)<f>-{z,x) = 6+(z,x)6-{z,x) 
and 

(5.52) 0±(A,x)0±(A,x) = <t>±(\x)0 ± {\,x), A€ct (2) . 



It remains to show that 0±(z, x) — 4>±(z, x). This is equivalent to showing that 

_ 0+(z,x) _ <p-(z,x) 
P[Z ' X) - 4> + {z,x) 9_(z,x)- 
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is equal to one. This function is well defined as long as <j>+{z, x) ^ 0. If 4>+(z, x) = 
this is still true as long as 9- (z, x) has no singularity (which is the case for z da- 
by Lemma [5.4p and does not vanish. But for z ^ da- the case <fi+( z , x) = 9- (z, x) — 
implies W(z) — 0, that is, z G ad U da+. Hence we will avoid such cases and 
suppose x ^ X, where 

X:= |J {x\$+{\,x)=0}. 

Fix x E M. \ X . Our aim is to show, that 

(5.53) p{z,x) = 1, z G C, x i X, 

By Lemma 13. II and Corollary [STB] the functions <j)±{z,x) and 9±{z,x) are holomor- 
phic in C\a± and hence p(z, x) is holomorphic on C\a^ with continuous limits up 
to the boundary away from da^ . By (|5.52[) the limits from different sides match 
up and so p(z, x) is even holomorphic on C\ da^ 2 \ Furthermore, arguing as before, 
one sees 



-(z,x) 9 + (z,x) 



<P+(z,x) 4> + (z,x) 



= O 



1 



E G da A 



that is, the singularities near E G da^ are removable and so p(z, x) is entire with 
respect to z for all x ^ X . Finally, p(z, x) — > 1 as z — > oo by item (iii) of Lemma l2.ll 
respectively Lemma T5. 41 In summary, f|5 . 53[) holds, that is, 

(5.54) 9±(z,x) = 4>±(z,x) 

for all x ^ X. But since the set X is discrete, this even holds for all x G R by 
continuity with respect to x. 

Finally, (|5.16|) shows that H+9 T (z, x) — z9 T (z, x), that is, q+(x) =q-(x). More- 
over, (|5.16|) and (|5.54|) imply, that 



T T (A)^±(A,x) = T (A,x) +i? T (A)0 T (A,x), 

and from (j5~50j) . (l5~54|) . ([3~16|) . and (|3~2T|) we conclude that ([3~TT]) is valid. Therefore 
the data S are the scattering data for the Schrodinger operator with the potential 
q(x) = <i-(x) — <7+(x). This finishes the proof of Theorem 15.31 □ 

Finally, observe that the second moment in condition IV from Lemma 14.21 can 
be replaced by any higher moment. In fact, introduce 

IV*. The function F±(x,y) is continuously differentiable with respect to both vari- 
ables and there exists real-valued function q±(x), id, with 

/>±oo 



± 



(1 + |x|")|g-|-(a:)|(ix < oo, for all a G 



such that 



\F±(x,y)\<C ± (x)Q± (x + y) 



Tx F±{ ^ v) 



<C ± (x) 



q± 



Q±{x + y) 



± 



dx 



F±(x,x) 



(1 + \x\ n )dx < oo, 
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where n — 2, 3, 4, ... , 



/•moo 

Q±(x):=±j x \q±(t)\dt, 



and C±(x) > is a continuous function, which decreases monotonically as x — » 
±oo. 

Then, proceeding literally as in Theorem 15. 3[ we obtain 

Theorem 5.8. Let the scattering data S, defined as in (|4.22|) . satisfy conditions 
I— III, and IV*. Then each of the GLM equations (14. ip has a unique solution 
K±(x,y). The functions q±{x) defined by (|5.7|) . coincide: 

<7-0) = q+( x ) = q( x ) 

and satisfy 

r±oc 

(5.55) ± \q(x) -p±(x)\(l + \x\ n )dx < oo, aeR. 



6. The Korteweg-de Vries equation with steplike finite-gap initial 

DATA 

In this final section we will outline in what way the inverse scattering transform 
method can be used to study the initial-value problem for the Kortcweg de Vries 
equation with steplike finite-gap initial data. Note that the Cauchy problem for 
steplike constant initial data is studied in [5J HU] ■ In the case of a constant back- 
ground on the right half-axis and periodic finite-gap background on the left one, 
this problem is partly considered in [T2] . The existence of the solution of the KdV 
equation with general finite-gap steplike potential as initial data seems to be an 
open problem and is not a subject of the present paper. Here we only discuss a 
possible approach. 

Let q(x) be a smooth function, satisfying condition (|5 . 55[) for n — 5 together 
with its derivatives: 

/■±oo 

(6.1) ± \q {k) (x)-p { ±\x)\(l + \x\ 5 )dx <oo, aeR, fc = 0,l,...,7. 



(6-2) ^7-6«- + — = 0, 



Here p±(x) are some finite-gap potentials. 
Consider the initial value problem 

du du d 3 u 

dt dx dx 

(6.3) u(x,0)=q(x), iGl. 

One can ask to find a unique smooth solution of this problem in the domain \t\ < T, 
satisfying conditions 

/■±oo 

(6.4) sup ± / {1 + \x\ 2 )\u{x 1 t) - u±{x,t)\dx < oo, 
\t\<T Jo 



/•±oo 

(6.5) sup ± / (1 + M) 

\t\<T JO 



d k u(x,t) d k u±(x,t) 



dx < oo, k = 1, 2, 3, 



dx k dx k 

where the functions u±(x, t) are the finite-gap solutions of equation (|6.2p with initial 
data u±(x, 0) = p±(x). We will proceed by a standard scheme 

u(x,0) ~» 5(0) -» S(t) u{x,t). 
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The Lax pair, associated with the KdV equation has the form 

(6.6) P(t) = -4^3 + 6u(x, i)— + 3u x (x, t), 

(6.7) H(t) = --^ + u(x,t). 

Equation (|6.2p is then equivalent to equation <9 t i/ = [H,P]. Let H±(t), P±(t) be 
Lax pairs, corresponding to our backgrounds. Following the scheme proposed in 
[S], we will obtain the time-dependent GLM equation: 

Lemma 6.1. Let ip± (z, x, t) be the Weyl solutions satisfying (see, e.g., [17) ) 

(6.8) H ± (t)il> = zij>, P±m = -i>, i> ± {z,0,Q) = l. 
Then the inverse scattering problem (the GLM) equation has the form 

/>±oo 

(6.9) K±{x,y,t)+F ± (x,y,t)± K±{x,s,t)F±(s,y,t)ds = 0, ±y > ±x, 

J X 

where 

(6.10) F ± (x,y,t)=I R ± (\,0)i> ± (\,x,t)i> ± [\,y,t)d P± (\,0) 



' ■ 



T T {\, 0)| 2 ^±(A, x, i)V>± (A, y, t)dp T {\, 0) 



+ 1] ik (°) 2 ^± ( A fc . *)^± y> *)• 

fe=l 



To prove that the problem (|6.2p - (|6.3p has a solution in the class (I6.ip - (I6.5[) 
one has to check, that time-dependent scattering data satisfy conditions of Theo- 
rem [231 Conditions I, II, III (a) can be verified directly from the time evolution 
for scattering data 



R±(\ t) = R±(X, 0) exp(a ± (X,t)-a ± (X,t)), 
T±(X,t) = T±(A,0)exp(a±(A,i) - a T (X,t)), 
7 fc ± W 2 =7fc(0) 2 exp(2a±(A fc ,i)), 

where (cf. [El US]) 

ct±(A,t)= jT (2(u ± (0, s) + 2A)m ± (A, s) - 9u± ^ s) \ ds 

and m±(X,t) are time-dependent Weyl functions, corresponding to background op- 
erators H±(t). Note that condition III, (b) means, that when E a v and lij(t) 
gets close to E, then R±(E,t) changes its sign. This effect was explained in [S]. 

To check condition IV one has to take into account the structure of the Weyl 
solutions ip±(X,x,t). As is known (see e.g., [I]), they admit a representation 
?/>±(A, x, t) — exp(±ifc±(A)a;)/±(A, x, t), where fc±(A) are the quasi-momentum maps. 
To estimate the parts Fd t ±(x,y,t) and Fh,±(x,y,t) of the kernel F± (cf. (|4.9p . 
(|4.15p , (|4.17p it is sufficient to use the Herglotz property of the quasi-momentum 
(Im(fc± (A)) > as A G M \ a±) and the fact, that the functions f±(X,x,t) are 
bounded as x € E, |t| < T and A ^ M±(0) U M±(0), A > inf{^,^, Ai} - 1. To 
estimate F rt ±{x,y,t) we integrate the first summand in (|6.10p twice by parts with 
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respect to the quasi-momentum variable k± and prove, that the boundary terms 
vanish. This approach fails only in points of the set [da_ U da^) D da^ (the 
points of type E$ in our example). Here one has to use the approach, developed in 
[20l Proposition 2.7]. This way one arrives at the estimates 



that justify condition (|6.4[> . 



9F±{x,y 1 t) 



Ox 



C(t) 
\x + yr 
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Appendix A. Properties of the transformation operators and 

ESTIMATES FOR THE GLM KERNEL 

In this appendix we derive and thoroughly investigate the integral equations for 
the kernels K±(x, y) of the transformation operators. In addition, we will obtain the 
necessary estimates for them and their derivatives. This will allow us to simplify the 
necessary and sufficient conditions on the functions F± (x, y) (in comparison with 
[14j ) and to solve the scattering problem in the prescribed class of perturbations 

O- 

Most of the results from this section are in essence known or follow as in the case 
of a constant background (see, e.g., [T31 [14] , and in the discrete case [HIH1[2T]). We 
included them here (with proofs) to make our presentation self-contained. 

Let ip±(z,x) be the background Weyl solution (|2.9p . Set 

(A.l) J ± (z, x, y) = ^ z ' M ~ ^ (Z ' V) 

and 

(A.2) q ± (x) = q(x)-p ± (x). 

Then the Jost solutions (|3 . 2[) satisfy the integral equation 

/>±oo 

(A. 3) (j>± {z,x)=i>±{z,x)- I J± {z, x, y)q± {y)4>± (z, y)dy. 

J X 

If we substitute formula (|3.2p into this equation, multiply with tp±(z, s)g±(z), and 
integrate over the set cr^.' 1 , using the inverse Fourier transform (|5.1ip . and taking 
into account that K±(x, y) = 0, ±x > ±y, we obtain 

/>±oo /> 

(A.4) K±(x,s)+ dyq ± {y)j> J± (A, x, y)ip± (A, y)4>± (A, s)dp± (A) 

/■ioo r ± oo p 

±/ dyq ± (y) dtK±(y,t)f J±(A, x, y)V±(A, t)^±(A, a)dp±(\) = 0. 

Set 

(A.5) r±(x,y,t,s) = T f ^±(Kx)^ ± (\,y)^ ± (\,t)^ ± (\,s)g ± (\)d P± (X), 
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where the integral has to be understood as a principal value. 

Then substituting (I2T4]) . (f2TT1) . (jA.lj) . and (|A.5|) into (|A.4|) we obtain 

/■zboo 

(A.6) K±(x,s) + / (r±(x,y,y,s) -T±(y,x,y,s))q±(y)dy 



p±OQ />±00 

±/ dyq±(y) K±(y,t) (T±(x,y,t,a) -T±(y,x,t, a)) dt = 0. 

J a: J y 

Consider now the function in (|A.5j) . From (|2. 13|) it follows that 



(A.7) r ± (x,i/ s t,s) = -r ± (i/,a!,«,t). 

Our plan is to evaluate the integral in (|A.5j) using the Jordan lemma. The only 
poles of the integrand in (|A.5[) are at the band edges and hence we introduce 

r± 

(A.8) f±(E,x,y) = Hm (j[{z - »f))^ ± {z,x)$ ± (z,y) 

Z ^ 3=1 

and 

2r± 

(A.9) D±(x,y,t,s)=±- ^ , P±(z) = [[{z-E j ). 

Note that D±(x,y,t 1 s) is a continuous and bounded function with respect to all 
variables. 

Now suppose ±(x — y + 1 — s) > and take a closed contour consisting of a large 
circular arc together with some parts wrapping around the spectrum <j± inside 
this arc at a small distance from the spectrum. Due to (|2.11[) . (|2.7p . and (hi) of 
Lemma [2TTT it follows that 

9±(z) a i/>±(z,x)#±(z,vW±(z,t)$±(z,a) = O (~) e ±V*(*-»+t-.) 

as z — > oo. In fact this holds on the entire circle since the neighborhood of the 
positive real axis can be handled as above. Hence one can apply Jordan's lemma to 
conclude that the contribution of the circle vanishes as its radius tends to infinity. 
Shrinking the loops the integral converges to 

(A. 10) T±(x,y,t,s) = D±(x,y,t,s), for ± [x - y + 1 - s) > 0. 

Note that f±(E, x, y) are real, and f±(E, x, y) — f± (E, y, x). Thus, the D±{x, y, t, s) 
is also real, 

(A.ll) D±(x,y,t,s) = D±(y,x,t,s). 

Now let ± (x - y + t- s) < 0, that is, ±(y - x + s - i) > 0. Then (|A.7j) . (|A.10j) . 
and (jA.lip imply 



T±(x,y,t,s) = -D±(x,y,t,s) = -D±(x,y,t,s), ±(x-y + t-s) < 0. 
Therefore, 

(A.12) T ± (x,y,t,s) =D±(x,y,t,s)sigD.{±(x-y + t- s)). 

Property (|A.11[) implies that the domain, where in the first integrand in (|A.6|) does 
not vanish, is 

(A. 13) sign(±(x - s)) = -sign(±(2y - x - s)), ±s > ±x. 
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In the second integral the domain of integration is 

(A. 14) sign(±(x -y + t-s)) = - sign(±(y -x + t- a)), 

with ± s > ±x, ±i > ±y > ±x. 

Solving (|A.13|) and (|A.14[) we arrive at the following result. 

Lemma A.l. The kernels K±(x,s) of the transformation operators satisfy the 
integral equation 

/>±oo 

K±(x, s) = -2 / q± {y)D ± (x,y,y, s)dy 

/£+£ 
/>±oo />s±y=pa: 

(A.15) t2 / dy D ± (x,y,t,s)K±(y,t)q±(y)dt, ±s > ±x, 

J x J 5±x=f y 

where D± are defined by (|A.9|) . 

Set s = x in (|A.15|) . Then the second summand vanishes, because we have our 
integration inside the domain ±t < ±y, where K±(y,t) = 0. Thus 

/>±oo 

(A.16) K±(x,x)=-2 q±(y)D±(x,y,y,x)dy. 

But, as is well-known (see, e.g, [13 Eq. (1.84)] or [22J Chapter 8]) 



r± r± 



(A.17) i> ± (z,y)i> ± (z,y) l[(z - M ±) = TJ(^ - M*(v))> 

j'=i i=i 

where /x^- (y) are the Dirichlet eigenvalues corresponding to the base point x = y 
(rather than x = 0). Combining (|A.5|) and (|A.10p we obtain 



D± (x, y,y,x) = D ± (x,x, y, y) 

(a, 8 , . ± I v ReS£ ns^fw) 



The integrand in (|A.18[) is meromorphic in C, thus, by the Cauchy theorem, we can 
compute the residue at infinity and obtain 

From (|A.16|) we conclude that 

(A.19) if±(a: s a:) - ±~ (q(t) - P±(t))dt. 



This formula justifies formula (|5.7[) under the condition, that the transformation 
operators kernels are differentiable. 

Lemma A. 2. Let 

/■±oo 

(A.20) Q±(x):=± \q±(t)\dt, q±(x) = q(x) - p±(x). 
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Then K±(x,y) has first order partial derivatives with respect to both variables. 
Moreover, for ±y > ±x the following estimates are valid 



(A.21) 
(A.22) 



\K ± (x,y)\<C ± (x)Q±(x + y), 



dK ± (x,y) 



Ox 



dK ± (x,y) 



dy 



<C±(x) 



q+ 



x + y 



+ Q±{x + y) 



where C±(x) are positive continuous functions for x G R which decrease as x — > ±00 
and depend on the corresponding background data and on the first moment of the 
perturbation. 

Proof. We restrict our considerations to the "+" case only and omit "+" in what 
follows. We will follow the scheme of the proof of [26l Lemmas 3.1.1, 3.1.2]. Intro- 
duce the following change of variables in (1A.15|) : 

(A.23) y + t=:2a, t-y=:2/3, x + s=:2u, s-x=:2v, 

then from (|A.15j) we obtain (see [26l Lemma 3.1.1]) 

/>oo 

(A.24) H{u,v) = -2 q(s)D 1 (u,v,s)ds 

J U 

poo fV 

-4 / da I q{a~ P)D 2 (u,v,a,f3)H{a,f3)dp, 

Ju JO 

where we put 

H(u, v) = K + (u — v, u + v), Di(u, v, s) — D + {u — v, s, s, u + v), 



(A.25) 



I?2(u, v, a, 0) = D + (u — v, a — f3, a + f3, u + v). 



Functions D\ and D2 are bounded uniformly with respect to all their variables. 
Put C — 2max{max„ i „ jS \Di\,maXu lV a ,/3 I-D2I} an d apply the method of successive 
approximations (see [26l Lemma 3.1.1]). We arrive at the estimate 



\H(u,v)\ < C(u~v) / \q{x)\dx 



(A.26) 

with 

(A.27) 



from which (|A.21jl follows. To obtain (| A.22|) . observe that the first partial deriva- 
tives of D\ and D2 exist (sec (|A.8|) . (|A.9|) and (|A.25j) ) and are bounded with respect 
to all variables. Thus, 



C(u) = Cexp C / Q(2t)dt), C > 0, 



dH{u,v) 
du 

= -4 



(A.28) 



— 2q(u)Di(u, v, u) = 

roc PIT) ( \ 

q{u-fi)D 2 (u,v,u 1 P)H{u 1 fi)dl3-2 / q(s) ° l^lll ds 

Ju ou 

*da [ V q(a-P) dD2{u > V ' a ' P) H(*,P)d(3, 
./n OU 
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dH(u, v) 



dv 



(A.29) 



9 — i(u, v, s ) ^ g _ 2 f q(a — v)D2(u,v,a,v)H(a,v) da 
dv J u 

da [ V q (a-/3) dD2{U > V ' a, ® H ^P)W 



The function Q(u) — J°° \q(x) \ dx is positive, monotonically decreasing, and satis- 

fies Q(-) e LHo.oo), B6R Since |A|, |^£H |^| <Ci, t = l,2, (1X261) applied 
to (fOg)) and ([X29]) implies 



dH(u,v) 



dv 



— 2q(u)Di(u, v, u) 



+ 


dH(u,v) 




dv 



< C(u-v)Q(2u), 



where the function C(u) is of the same type as (|A.27j) . with a different positive 
constant Ci depending on the background data. From this, (|A.23j) . and (|A.25jl . the 
estimate (|A.22[) follows. □ 

With the help of this lemma we can now derive several estimates for the GLM 
equation. 

Lemma A. 3. The kernel F± (x, y) of the GLM equation (|4. 1[) has first order deriva- 
tives with respect to each variable. Furthermore, for ±y > ±x it satisfies 



(A.30) 



\F±(x,y)\<C ± (x)Q ± (x + y), 



(A.31) 



dF±(x,y) 



dx 



dF ± (x,y) 



dy 



<C±(x) 



q± 



x + y 



+ Q±(x + y)\, 



where the functions q±{x) and Q±(x) are defined in (|A.20j) . Here C±{x) are positive 
continuous functions which decrease as x — > ±oo. Moreover, 



(A.32) 



± 



±00 



dF± (x, x) 



dx 



< oo, Va € 



Proof. Again we restrict our considerations to the "+" case only and omit "+" in 
what follows. Set Qi(u) = Q(t)dt. Due to condition (|1.3p . the functions Q{x) 
and Q\(x) satisfy 



QiO)di<oo, / Q(t){l + \t\)dt < 



(A.33) 



Observe also, that the kernel F(x,y) of the GLM equation (|4.1[) is symmetric: 
F(x,y) = F(y,x). From (jlTTj) and (|A.21|) we see, that 



\F(x,y)\<C(x)(Q(x + y)+ / Q(x + t)\F(t, y)\dt 



Since Q\{x + t) > Qi(2x), and C(t) < C{x) as x < t, then Gronwall's inequality 
implies (|A.30p with 



C{x) = C 1 C(x)exp(CiC(x)Qi(2x)), d > 0. 
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Differentiating (14. ip with respect to x and y implies 

/>oo 

(A.34) \F x (x,y)\<\K x (x,y)\ + \K(x,x)F(x,y)\ + \K x (x,t)F(t,y)\dt, 

J X 

POO 

(A.35) F v (x,y) + K y (x,y) + / K(x,t)F y (t,y)dt = 0. 



The functions Q{x), Qi(x), C(x), C(x) are monotonously decreasing and positive. 
Furthermore, 



q± 



X + 1 



+ Q{x + 1) Q(t + y)dt < (Q(2x) + Q 1 {2x))Q{x + y), 



and hence the estimate (|A.31[) for F x follows (with some other positive continuous 
decreasing function C{x)) from (|A.30fl . (|A.2ip . (|A.22jl and (|A.34|1 . The same esti- 
mate for F y can be obtained from (|A.35I) and Lemma IA.2I by using the method of 
successive approximations. 

It remains to prove (|A.32[) . To this end consider (|4. 1[) for y — x and differentiate 
it with respect to x: 



dF(x,x) dK(x,x) 



-K(x, x)F(x, x)+ / (K x (x, t)F(t, x) + K(x, t)F y (t, x)) dt = 0. 



dx ' dx " " 1 ' ' 
Formula (|A.19jl implies Next, by (|A.30j) and (|A.21j) . we have 

\K{x,x)F(x,x)\ < C(a)C(a)Q 2 (2x) for x > a, 
where / a °°(l + x 2 )Q 2 (2x)dx < oo. Moreover, by (|A~3T|l and ((A~22| . 

\K' x (x,t)F(t,x)\ + \K(x,t)F y (t,x)\ <4C(a)C(a){\q(^-)\Q(x + t) + Q 2 (x + t)}, 
and together with the estimates 
dxx 2 I Q 2 (x + t)dt < 



\x\Q(2x)dx sup / \x + t\Q(x + t)dt < oo, 

x>a J x 



(^±i)\ Q (x+t)dt 



< 



< 



Q(2x)dx sup 

x>a 



(^-)\(l + (x + t) 2 )dt 



< OO 



we arrive at (|A.32[) . 



□ 



Remark A. 4. Note that the results of this lemma are in some sense invertiblc. 
Namely, if we start with properties (I4.18|) - (|4.21|) of F±(x,y), then, using (|4.ip and 
the same considerations as in Lemma [A. 3| we obtain (|5.1|) . (|A.22|) . and (15.21) . 

Appendix B. Proof of Lemma [B~4l 



We introduce the local parameter r = \J z — E in a small vicinity of each point 
E e da± and set y(z, x) = -g^y{z 1 x). Since 4| (E) = 0, for every solution y{z, x) of 
the equation (|2.1[) . its derivative y(E,x) is also a solution of (|2.ip . In particular, 
the Wronskian W(y(E), y{E)) is independent of x. 
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For each x £ R in a small neighborhood of a point E € da± introduce the 
function 



(B.i) 4>±A*>x) 



ip±(z,x), E e da ± \ M±, 
ti/)±(z,x), E £ M ± . 

Lemma B.l. Let the function 'ip± l E{z, x) be defined by formula (|B.1|) . Then 

: lim r—r, 

z^E 2g ± {z) 

where a = — 1 if E £ da± \ M± and a — 1 if E £ M± . 

Proof. We begin by recalling (see, e.g., [TTl Eq. (1.73)]) that the Weyl m-functions 
can be written as 



(B.2) w($±,eW,-^$± iB (E)) =±Um 



(B.3) m±(z) = — ^ , 

where G±(z) = \F' ± (z,0) and F±{z) = F±(z,0) with 

2r± r± 

(B.4) P±(z) = *[[{z-Ef), F ± (z,x) = l[(z-nf(x)). 

3=0 j=l 



Furthermore, observe that Z f\/P± (z) — rf±(z), where f±(z) is holomorphic near 
z = E. 

Then in the case E £ da±\M± (where ^±,s(A, x) = tp±(X 7 x)) we have F±(E) ^ 
0. Then, from l[2~3]) . we have i/)±(E,0) = 1, ^±(-E,0) = m±(E), ip±(E,0) = 0, 



^,0) = m ± (^ = ^| = T hm^, 

r±(E) z^-E tL±(z) 

and the first claim follows. 

In the second case we have £7 € M± (where $±,e(A, x) — Tip±(A,x)) we have 
F±(E) = and hence F±(z) = t 2 F ± (z), G±{z) = r 2 G±(z), where F±(z) and 
G±(z) are holomorphic near z — E with F±{E) ^ and G±{E) ^ 0. Hence we 
have 

i>±AE,o) = o, ^ (£?,o) = -^§, 

and the second claim follows. □ 

From (|2.5[) . (|2.10[) and (|2.9p we observe the following property of the Floquet- 
Weyl solutions 



Remark B.2. Let E £ M± and x £ R. Then ij)± tE (E,x) = -i>±,E( E '> x ) if -E is a 
left band edge from a± and ip± t E{E, x) — ip±, E (E, x) if E is a right band edge. 

It is straightforward to check that this property is also inherited by the Jost 
solutions. Again we abbreviate 



(B.5) <£±,e(A,z) 



4>±{\x), Eeda±\M±, 
T(f>±(X,x), E £ M±. 
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Lemma B.3. We have 4>±(E, x) = <f>± (E, x) for E G da ± \AI±. If E e M±, then 
by (|B.5[) 4>±,e{E,x) = —<j>± t E(E,x) when E is a left edge of a band from a± and 
4 ) ±,e{E 1 x) = 4>±,e(E, x) when E is a right edge of a band. 

Lemma B.4. The function W(z) is continuous on the set C\a up to the boundary 
er u U a 1 . It can have zeros on the set da U (da]P H da^) and does not vanish 
at the other points of the set a. If W(E) — as E G da U (da^ fl da_), then 
W(z) = Vz~E{C{E) + 0(1)), C{E) £ 0. 

Proof. Continuity of W up to the boundary follows from the corresponding property 
of <j)±(z,x). We begin with the investigation of the possible zeros. 

Let A G int(cr( 2 )) := a^ \da^ and suppose W(X ) = 0. Then c/) + (X Q ,x) = 
cJ>-( X ,x) <fi+( X ,x) = c4>-(X ,x), i.e. W(4>+,4^) = \c\ 2 W(4>-, JZ). But by 
(|2.1ip and (|2. 14[) we have signg + (Ao) = — signg_(A ), contradicting (|2.6|) . 

Let Ao G 'mt(a^) and W(Xq) = 0. The point Ao can coincide with a pole 
/i G A/zp. But ^>±(Ao,a;) and 4>±(Xq,x) are linearly independent and bounded, and 
4> T (x, A ) € K (see (|3~7)) ). If W{X ) = 0, then T = cf <f>± = 0±, that implies 
W(0±,0±)(Ao) = 0, which is impossible. 

In the general mutual location of the background spectra one can meet the 
case when Ao = E € <9ct ( - 2 - ) n int(er±) C mt(<r) (that is a point like E5 in our 
example). If W(E) = 0, then W(<f>±, 4> t ,e){E) = 0, where (j) T ,E is defined by 
(IB. 51) . But according to Lemma IB. 31 the values of cj)^ t E(E, ■ ) are either pure real 
or pure imaginary, therefore W(4>±,4 > t,e)(E) = 0, that is, (j>±(E,x) and <j>±(E,x) 
are linearly dependent, which is impossible at inner points of the set a±\ 

In summary, W(X) ^ for A S int(cr) \ (da [ l ] n da { + ] ) which finishes the first 
part and it remains to investigate the order of zeros. 

Let E G da U (da^ D da_) (these are the points of type E\, E3 and E± from 
our example). The function W(A) is continuously differentiable with respect to the 
local parameter r. Since at E we have -4^(5+8_)(E) = 0, the function W(<j) + ,(j)J) 
has the same order of zero at E as W(X). But if 5(E) ^ 0, then 4j-5±(E) = and 
if L(E) = 5+(E) = 0, then ^(t~ 2 5+ L)(E) = 0. Therefore jpW(E) = if and 
only if f W(4> +i e,$-,e) =0/ 

To simplify notations, we will just write <j>± := 4>±^e until the end of this proof. 
Again we have consider all possible cases for the mutual location of the Dirichlct 
eigenvalues. 

First let E G da^ n da (a point of type E2 in our example) and let E 
(M+ U M_). Then W(E) = if and only if (see l|5^]) ) = W^(0+, <M = 0, 

that is, <fi±(E, ■) = c±<j) T (E, •), c_c+ = 1, c_,c+ G R. The derivative of the 
Jost solution with respect to r is again a solution of equation (|3.ip . Therefore, by 
Lemma IB. 11 

(B.6) W(E) = W(h,^-) - W(j>-,<l>+) = c_ W(<p+^+) - c+ W{<j>-,<l>-) 
= c_ W{i>+, V+) - c + W(^_, V-) = -(c+ d_ + c_ d+), 

where 

d+ = lim ■ 

a^b 2 ff± (A)v / A^TB 
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We see from jUBJ that d± G iM + \ {0} if E is a left edge of a and d± G R+ \ {0} 
if E is a right edge of a. Since signc_ = signc+, this finishes the case E G 
da^ nda\(M+UM-). 

The same arguments are valid for E G da^ (1 da (1 M + n M_ and PF(E') = 0. 
Then <j)±(E, ■ ) = c± (f> T (E, ■ ) and using Lemma [B~3l we conclude that c_, c + G M, 
sign c_ = sign c + . By Lemma IB. II 



±W(hJ-)(E) = c + lim ^EA + c _ & ^ o. 

The cases, when one of the Jost solutions is bounded in the edge of spectrum, 
and the another one is unbounded, are more subtle. For example, let W(E) = for 
E G (M + ndcr (2) n<9er)\ (M- HM+) and let £ be a right band edge (like the point E 2 
of or example when /i+ = E% and /i_ ^ E%)- Then by Lemma lB.3l c!)+ fx. E)gl and 
4>-(x,E) = (/)_ (x, E) G M. Also ^ if and only if £W($+, (j>-)(E) ^ 0. 

Therefore we have = c + <j!>_, <t>— = c -4>+ an d signc+ = signc_. Thus, by 
Lemma IH31 j-M^, (j>-)(E) = C-d+ — c + eL, where 



d+ = lim ^ A /x f eL\ {0}, d_ = lim 1 — S M+ \ {0}. 

+ a^b 2i 5+ (A) ^ l ; ' A^BiyX^B2i 5 _(A) 

The case, when E G (M_ H <9cr( 2 ) n <9ct) \ (M+ n M_) and £ is a right band edge is 
analogous. 

Now, let E G (M + n 9cr (2) n da) \ M_ and let J5 be a left band edge. Then by 
Lemma lB.31 d>+ G iM and c+,c_ G iK, c + c_ = 1, that is, sign(ic + ) = — sign(ic_). 
Furthermore, sign y/X — E > since X > E, and signg + (A) = signg_(A). There- 
fore, 



0-)(S) = ic_ lim — —J - i c + lim 1 — ^ 0. 

Unlike the case of one and the same background p+ (x) = p_(x), where W(X) ^ 
for A G int(cr), we could have W(X) = for A G int(cr) in our steplike situation. 
The points under consideration are points of the set da^ n da^ . This case can 
be treated in the same way as the case E G da {2) n da. Namely, if E $ (M_ U M+) 
then observe that one of the two summands in (|B.6|) is real and the other one 
is imaginary. The same is valid for the case E G (M- n M+). Now let S G 
da^ nda^ nM + \(M + flM_) and let £ be a right band edge of ct^ (i.e. a left band 
edge of a^). Then by Lemma MM $+(x, E) G iM and $-(x,E) = <p-(x,E) G M. 
Moreover, VF = if and only if W((/)+,(j)-) = with the same order of zero. 
Therefore, c+, c_ G iR and by pO) and 

^W(<£+,<M = c_VF(^ + ,^+) -c+^(4->-) 

= — fic_ lim h ic+ lim ") 7^ 

27T V a^_e g + x^E rg- J 

because the first summand is imaginary and the second one is real. All other 
combinations can be treated similarly. 

Finally, consider the case E G da±~* n da. Let, for example, E G da^l n da 
(the point E4 in our example). Since E G M \ cr_ in this case, we have 6 -(E) 7^ 
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and one has to study zero of the Wronskian W(4>+,<j)-). But by (|2.8j) and (|3.2j) 
W§-, = W$-, £pi>-) = 0. Therefore, 

^W($ + ,$-)=c-W(Jp$ + ,$+) ^0 
by Lemma [B. II □ 
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